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Abstract 

The wave function and binding energy for shallow donors in GaAs are calcu- 
lated within the tight binding (TB) approach, for supercells containing up to 
two million atoms. The resulting solutions, coupled with a scaling law, allow 
extrapolation to the bulk limit. A sharp shallow-deep transition is obtained 
as the impurity perturbation increases. The model allows investigating the 
quantitative consistency between the effective mass theory and the TB for- 
malism. Although the calculated binding energies are in excellent agreement, 
anisotropies and the overall decay obtained in the TB envelope function can 

not be afforded by the hydrogenlike effective mass prediction. 
PACS numbers: 71.15.Ap, 71.55.-i, 71.55.Eq 
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Electronic states in condensed matter systems may exist in either localized atomic-like 
states or in delocalized running- wave states [|1|. Such dualistic behavior has led naturally 
to two well-established frameworks in which to describe the electronic behavior of solids: 
Real space pictures privilege the atomic aspects of electronic behavior, whereas Bloch states 
provide a foundation to describe the wave-like properties. For an ideal periodic crystal, 
these two views were reconciled several decades ago, when Wannier showed that a set 
of orthogonal "atomic" functions localized at the crystal's atomic sites may be obtained 
from the Bloch eigenfunctions, and that either may be taken as the basis set for the general 
description of electrons in solids. Any electronic theory which adopts a basis set localized 
in real space is generally referred to as following the tight-binding (TB) approach. In some 
practical cases, it is also possible to utilize a basis set localized in reciprocal space about an 
appropriate point in the Brillouin zone. Such an approach then gives rise to the so-called 
Effective Mass Theory (EMT) 0. EMT is the most successful comprehensive approach for 
the description of electronic excitations in insulators, of transport and optical properties in 
semiconductors, of dynamics of electrons in metals [Q. 

A crowning achievement, dating as far back as the 1950's, of the EMT is undoubtedly 
the quantitative understanding obtained of shallow impurity states in semiconductors |5|-0 . 
A complementary TB theory for deep impurities was presented later by Hjalmarson et al, [§. 
The choice of approach in each case is defined naturally, given the extended nature of the 
shallow levels described within EMT versus the localized nature of the deep states considered 
in the TB scheme. A basic limitation of the EMT, however, is that the approximation 
(i.e., a single dominant k-point) breaks down as the impurity state becomes increasingly 
delocalized throughout the Brillouin zone, such as the case for deep levels. Similarly, a 
practical limitation of the TB approach is that until now it could not be extended to the 
shallow impurity limit due to the slow decay of the defect wave function which requires 
extremely large supercells for meaningful calculations. As a result, for decades a gap has 
existed between two of the most widely used theories for impurities in semiconductors: The 
TB approach was confined to deep states, whereas the EMT was used for shallow states. 
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In the present work we demonstrate that the gap between shallow and deep levels may 
be bridged via the TB formalism in conjunction with extrapolation schemes for the binding 



energies |T2| and for the envelope functions. Our proposed atomistic description provides 
a flexible way to explore and test key predictions of EMT for shallow levels. In particular, 
the following points are addressed and explained: (a) The values of ionization energies, (b) 
the independence of these ionization energies on the impurity species, and (c) the shape of 
the envelope function for the defect state. We show that wave functions may substantially 
deviate from the hydrogen-like predictions of EMT. Besides the obvious academic interest of 
unifying shallow and deep impurity level descriptions within a single theoretical approach, 
problems of current interest in nanostructured materials demand an atomistic description for 
which the TB approach coupled to the theoretical tools adopted herein may provide a reliable 
alternative to EMT and to ab-initio treatments based in the local density approximation. 

In the EMT, it is assumed that the impurity wave function is highly delocalized in real 
space, which in turn implies a strong localization in k-spa.ce. The electronic properties of 
the host material are thus described by a few parameters related to the energy bands near 
the set of A;-points in consideration. In particular, effective masses are defined from the 
bands' curvatures at these points. The simplest case is that of a nondegenerate conduction 
band with an extremum at F (k = 0), corresponding to the Bloch eigenstate ipc{r) and 
eigenvalue ^^^(r). The cubic symmetry of the lattice implies that the nondegenerate band 
is isotropic around this extremum, leading to a single effective mass m*: Ec(k ^ 0) = 
Ec{T) + /i^k^/2m*. Following the original EMT formulation 0, the perturbation potential 
for donor impurities is 

U{r) = -e^er, (1) 

where e is the static dielectric constant of the host. Within a set of additional approximations 
0, the impurity wave function is written as \E'(r) = F{r)ipQ{r), and the eigenvalue problem 
for the one-electron Hamiltonian H = Hq + U{r), where Hq describes the perfect host 
material, maps into an hydrogenic equation 



[-h' v' /2m* + f/(r)]F(r) = [E - Ec{T)]F{r) (2) 



leading to the ground state envelope function F{r) = (l/vTra*'^) exp(— r/a*), where a* = 
ao£:(mo/m*) is the effective Bohr radius, ao=0.53 A is the Bohr radius of the hydrogen atom 
and rriQ is the free electron mass. The EMT basic assumptions are fulfilled if a* is much 
larger than the host lattice parameter. In this case, the donor binding energy for a single 
impurity is = Eh [m* /{e'^mo)], where i?/f=13.6 eV is the hydrogen ionization energy. 

For donor levels in semiconductors with a conduction-band minimum at F, the above 
model combines simplicity and physical intuition and moreover leads to a good agreement 
with experimental results. In the case of GaAs, e = 12.56 and m* /itlq = 0.068, leading 
to a* = 97.7 A (which is indeed much larger than the host lattice parameter acaAs = 
5.65A), and to E^ = 5.86 meV, which is in remarkable agreement with the experimental 



data for donor binding energies in this material |jT3[. It is intriguing that most donor levels 
in GaAs are essentially independent of the substitutional species, e.g. 5.91, 5.84 and 5.88 
meV respectively for C, Si and Ge replacing Ga, while acceptor levels vary by almost a factor 
of two for these same species substituting As [|l^ . We choose this extreme situation of very 



shallow donor levels, for which EMT predictions for the binding energy work remarkably 
well, to explore the limits and capabilities of the TB formalism within a recently proposed 
scheme |T2[ . 

We consider a single impurity placed in a large cubic supercell containing = 
8(L/aGaAs)^ atoms arranged in the zincblende structure, where L is the supercell edge length. 
Periodic boundary conditions are imposed, and supercells containing up to 2 million atoms 
(L = 64 X acaAs = 361.6 A) are considered. We use the sp^s* basis for the TB description of 
the electronic structure, with the bulk GaAs Hamiltonian taken from Boykin's parametriza- 



tion which includes first and second neighbor interactions and leads to an accurate 
value for the conduction band effective mass at F. Due to the s-character of the GaAs 
conduction band edge, spin-orbit corrections may be neglected in our calculations. The TB 



Hamiltonian is written as [12 
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H =Y1 hfjcl^Cj^ + {r,;)c\^Ci^ (3) 

where i and j denote the sites in the zincblende structure, /i and v denote the atomic orbitals 
and Ti is the distance from site i to the impurity site. The h^^ define all the on-site energies 
and first and second neighbor hoppings for the bulk material. The perturbation potential 
f/(rj) is described by Eq.([|), except at the impurity site (rj = 0), where it is assigned the 
value f/o, describing central cell effects according to the substitutional species. In the present 
calculations, f/o is kept as a free parameter. Estimates for this parameter are of the order of 
one to a few eV Note that the impurity potential in incorporates both the long- 

range component of the original EMT formulation 0] and the short- range perturbation {JJq 
restricted to the impurity site) of the original TB formulation p. 

The exact ground state wave function and binding energy for a donor level within a 
supercell of size L was obtained by minimizing the expectation value of 
where Sref is a reference energy chosen below the conduction band minimum with the 
TB wave function expansion coefficients {Cj^} taken as independent variational parameters. 
This procedure is equivalent to the exact diagonalization of H for the eigenvalue and eigen- 
function closest to Sref- The calculated binding energies as a function of f/o are given in 
Fig. |l|(a). It is clear that there is a range of the values of f/o for which El{Uo) is constant, 
indicating that the binding energy does not depend on the impurity species, in agreement 
with the EMT and with experiments A well defined shallow-to-deep transition occurs 
for f/o ~ 1.8 eV, above which the binding energy increases approximately linearly with 
Uq. Aiming at a better characterization of this transition, we calculate the orbital averaged 
spectral weight of the donor state, 

(k) = ^ E e^'-^^^~^^^a,Q„ (4) 

where Rj and Rj denote the position vectors of the i and j atomic sites. The calculated 
W (k) for k = r, X and L points in the fee Brillouin zone are given in Fig. |l](b). We visualize 
in this Figure the /c-space counterpart of the shallow-to-deep transition around f/o = 1.8 eV. 
For f/o ^ 1.8 eV, the wave function has an almost pure F character, since Vr(r) is very 
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{H - ErefY 



close to 1. This extreme localization in /c-space is consistent with the EMT assumptions, 
and indicates a highly delocalized wave function in real space. Increasing Uq beyond 1.8 eV, 
the dominant F character breaks down and the impurity state must be described in terms 
of the several k-point components, indicating delocalization in A;-space, localization in real 
space. This is consistent with the deeper character of the state obtained in Fig. |I|(a) in this 
range of f/o. 

In their original TB study, Hjalmarson et al noted that shallow levels are unbound 
by the central cell potential alone, and thus could not be obtained there. Our calculations 
are in complete agreement with this interpretation: We find that, in the whole range of 
Uq for which El{Uq) is constant, no bound state is obtained if the Coulomb potential is 
not included in the impurity perturbation. This result confirms a simple physical criterion 
characterizing shallow levels, and the reason why they do not depend on the substitutional 
impurity species, i.e. on f/o, since the binding is entirely due to the Coulomb attraction part 
of the perturbation potential. On the other hand, when deeper levels are obtained, f/o alone 
also leads to a bound state, and therefore the binding energy is sensitive to the impurity 
species. 

In the following we focus on shallow donors. Given that the binding energy is independent 
of f/o in this regime, we arbitrarily take f/o = 1.0 eV in the model calculations. This 
is a reasonable choice given the value of the Coulombic potential at the nearest-neighbor 
distance. Fig. |^(a) gives the donor ionization energy calculated at increasing supercell sizes. 
Results for presented there, even for the largest sizes, are still decreasing with L. This 
difficulty was overcome in Ref. |jl2| by assuming that, asymptotically, the convergence of 



to the isolated-impurity limit (L oo) should be exponential, namely 

EL = E^ + Ee~'^'\ (5) 

where E, A and the donor ionization energy for infinite cells E^o are taken as adjustable 
parameters. The validity of Eq.(^ implies a linear dependence of \n.[{Ei — Eoo)/E^] on 
the system size L. These are the circles in Fig. §(b). As indicated by the straight solid 



line, a linear behavior is obtained for L/a* ^ 1.5 with the extrapolated bulk donor binding 
energy E^o = 6.7 meV and A = 1.25a*. We note a much slower convergence of the donor 



shallow level with L as compared to previous calculations for acceptors ||T2[. Similar fits were 
obtained for Eoo = 6.7 ± 0.7 meV, while values of E^o outside this range lead to systematic 
deviations of the calculated data points from a linear behavior. The extrapolated binding 
energy E^o is indicated in Fig. §(a). Although the EMT value is slightly off the lower limit 
estimated numerically, our results do not indicate any significant discrepancy between the 



binding energies determined by the two approaches [|1^ . 

We now focus on the analysis of the donor wave function. Fig. |^ gives the TB envelope 
function squared (expansion coefficients squared, summed over the five orbitals for the cation 
sublattice), calculated for L = 48aGaAs ~ 

270A, function of distance from the impurity 
along the [100], [110] and [111] directions. At long distances, the wave function shows no 
angular dependence, as predicted by EMT. However, in the vicinity of the impurity there are 
noticeable anisotropies, which are highlighted in the inset. This is an effect of the crystalline 
environment that is automatically captured in an atomistic method such as TB. A more 
important discrepancy with respect to EMT is that the donor-state wave functions cannot 
be described by a simple exponentially decaying behavior. In fact, the line in Fig. |^ is a sum 
of exponential functions centered at each supercell-periodic replica of the impurity, decaying 
with the EMT Bohr radius of 97.7 A. The agreement is not good, and it cannot be improved 
by varying the Bohr radius. 

It is not surprising that EMT works well for the binding energies but not so well for wave 
functions. This can be seen as a manifestation of the variational principle: Small variations 
in the ground-state wave functions do not affect the ground-state energies in linear order. 
We test this hypothesis through a set of variational calculations in which trial wave func- 
tions are constructed by rescaling the TB expansion coefficients of the true wave function 
in such a way that they are constrained to an envelope composed of summing exponen- 
tial functions centered at each periodic replica of the impurity. In other words, in these 
test calculations, the TB wave function expansion coefficients {Cj^} are not independent 



variational parameters, but follow an exponential-decay behavior, corrected for the periodic 
boundary conditions. The decay length is now the sole variational parameter and it is then 
varied to minimize energy. We name these calculations "hydrogenic tight-binding" (HTB). 
The binding energies resulting from these constrained variational calculations also scale with 
supercell size as in Eq. (^, as shown by the triangles in Fig. |^(b). The extrapolated value, 
= 6.4 ±0.6 meV, is for all practical purposes indistinguishable from the full TB value. 
We also calculate, for different supercell sizes L, the optimum decay length ql, as illustrated 
in Fig. |. We find that the calculated lengths also converge exponentially: 

«L = Ooo + S e~^^'^ . (6) 

The extrapolated Bohr radius is Ooo = (93±5)A, entirely consistent with the EMT prediction! 
Therefore, even though the true wave function is not exponentially decaying, if such a 
constraint is imposed (as in EMT), the optimum Bohr radius agrees with EMT predictions 
and the binding energy is essentially correct. 

In conclusion, we study the problem of shallow donors in GaAs within the TB formal- 
ism. The use of large supercells (up to L = 361.6 A), together with a finite-size scaling 
procedure, allows an unprecedented test of EMT for shallow levels in semiconductors. We 
find good agreement between TB and EMT in the determination of the binding energy, but 
hydrogenlike wave functions do not provide an accurate description of the impurity state. 
In addition, by varying the central cell corrections to the potential, we confirm the EMT 
predictions that the binding energy is independent on atomic impurity species for shallow 
levels. Our results are inconsistent with the commonly accepted idea that it is straightfor- 
ward to describe bound states within atomistic approaches if the supercell is large compared 
to the Bohr radius (or typical localization lengths). For the binding energy, the value ob- 
tained for a cell at least 3 times the EMT Bohr radius is about 30% above the converged 
value. An extrapolation scheme is required to reach reliable results. Description of the wave 
function is even more subtle, due to the adopted boundary conditions. Periodic boundary 
conditions requires careful consideration of contributions from the impurity periodic images. 
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We demonstrate these inconsistencies and propose general and original ways to overcome 
them, applicable to any atomistic approach. 

We thank H. Chacham and F.J. Ribeiro for helpful discussions, and CNPq, CAPES, 
FAFERJ, MCT-FRONEX, and FUJB for financial support. 



9 



REFERENCES 
[1] P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] G.H. Wannier, Phys. Rev. 52, 191 (1937); Rev. Mod. Phys. 34, 645 (1962). 
[3] J.C. Slater, Phys. Rev. 76, 1592 (1949). 

[4] See for example J.M. Ziman, Principles of the Theory of Solids, Second Edition (Cam- 
bridge University Press, Cambridge, 1972); P.Y.Yu and M. Cardona, Fundamentals of 
Semiconductors, Second Edition ( Springer- Verlag Berlin, 1999), and references therein. 

[5] C. Kittel and A. H. Mitchell, Phys. Rev. 96, 1488 (1954). 

[6] J. M. Luttinger and W. Kohn, Phys. Rev. 97, 869 (1955). 

[7] S.T. Pantehdes, Rev. Mod. Phys. 50 797 (1978). 

[8] H.P. Hjalmarson et al, Phys. Rev. Lett. 44 810 (1980). 

[9] S.T. Pantehdes and C.T. Sah, Phys. Rev. B 10, 621; 10, 638 (1974). 
[10] T. B. Boykin, Phys. Rev. B 56, 9613 (1997). 
[11] J. G. Menchero and T. B. Boykin, Phys. Rev. B 59, 8137 (1999). 

[12] J. G. Menchero, R. B. Capaz, B. Koiller and H. Chacham, Phys. Rev. B 59, 2722 (1999). 

[13] O. Madelung, Editor, Semiconductors - Group IV Elements and III-V Compounds, Data 
in Science and Technology, Springer- Verlag (1991) 

[14] R. B. Capaz, G. C. de Araujo, B. Koiller, and J. P. von der Weid, J,. Appl. Phys. 74, 
5531 (1993). 

[15] Note the large energy scales involved in the problem, e.g. in Figs, ^(a) and |^(a), from 
which differences of tens of a meV may well be due to some degree of arbitraryness in 
the fitting procedures or small inaccuracies in our numerical variational scheme. 



10 



FIGURES 



> 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

1 

0.1 
1E-3 4 
1E-5 
1E-7 
1E-9 



-i — I — I — I — I — I — I — 1 — I — I — I — I — I — I — I — I — r 
o L= 339 A 



(a) 



o L = 226 A 
V L=113A 
A L = 90.4 A 



o 
o <f 



o 



-l — I — I — I — I — I — I — 1 — I — I — I — I — I — I — I — I — l- 




0. 



1 — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — 
8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 



U, (eV) 



FIG. 1. (a) Binding energy of the impurity level as a function of the central cell perturbation 
potential Uq for different supercell sizes L. (b) Calculated spectral weights for the F, X and L 
/c-points as a function of Uq- The solid lines are guides for the eye. Both (a) and (b) revel a clear 
shallow-to-deep level transition around C/q = 1-8 eV. 
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FIG. 2. (a) Tight-binding results for the binding energy of a single donor substitutional impu- 
rity in GaAs as function of supercell size L, with a central cell correction Uq = 1.0 eV. The dotted 
line represents the EMT value {E^ = 5.86 meV), which is in good agreement with experimental 
data. The circle at the lower left gives our extrapolated value for L — > oo, = 6.7±0.7 meV. (b) 
Exponential convergence of the donor binding energy for both TB and HTB calculations described 
in the text. The line is an exponential fit to the TB results (see text) with E^o = 6.7 meV. 
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FIG. 3. TB envelope function squared along the indicated lattice directions. The solid line is a 
sum of exponential functions centered at each periodic replica of the impurity, decaying with the 
EMT Bohr radius of 97.7 A. The inset illustrates the anisotropy of the TB wave function around 
the impurity site. 
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FIG. 4. Exponential convergence of the Bohr radius with the supercell size L in the HTB model. 
The dotted line corresponds to Uoo = 93 A. 
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